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1 Introduction 

The electromagnetic and strong interactions are well recognized as gauge theories. Nevertheless, 
compared to their elegant mathematical formulation, gauge symmetries play an perplexing role 
when we try to understand the angular momentum structure of the photon field and gluon 
field. On one hand, the spin and orbital angular momentum of light modes have been measured 
experimentally [1-3], and the gluon spin contributions to the nucleon spin are also believed to be 
measurable [4-6]. On the other hand, it is believed that a gauge invariant operator decomposition 
of the gluon spin and orbital angular momentum is not feasible, compared to a well established 
gauge invariant decomposition of the quark spin and orbital angular momentum [7]. 

Recently, Chen et al. [8] proposed that the gluon spin and orbital angular momentum can be 
identified gauge-independently by decomposing the gauge field into its physical part ^phys 
its pure gauge part ^pu^e- The decomposition of Chen et al. shows the appealing feature that 
the gauge covariant derivative operator in their construction appears to depend solely on the 
pure gauge field, so it can be regarded as the closest akin to the interaction-free canonical ones. 
Moreover, Chen et al. argue that the Coulomb gauge construction plays an especially superior 
role, mainly because the Abelian gauge field can be uniquely decomposed into its transverse 
part and longitudinal one. The proposal of Chen et al. has inspired lots of considerations on 
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the feasibility of gauge invariant operator descriptions of the gluon spin [9-13]. However, there 
are also some concerns about the decomposition of Chen et al., such as the frame- dependence 
problem in the Coulomb gauge case and the light-cone gauge case [14], the uniqueness problem 
of identifications of the gluon spin operator [15, 16], whether the requirement of gauge invariant 
operators is necessary [17], and whether the decomposition of Chen et al. can provide new un- 
derstandings on the decomposition of the nucleon spin beyond the framework of gauge invariant 
extensions (GIEs) [14]. In this paper, we discuss the feasibility of gauge invariant descriptions 
of the gluon spin and attempt to shed light on these concerns from several different aspects. 

Firstly, we propose that gauge invariant expressions of angular momentum current can be 
constructed from summations of two kinds of Noether currents: one kind of current is induced 
by the Lorentz transformation and the other one is induced by the gauge transformation. These 
results are presented in section 2. Secondly, in section 3 we propose novel expressions for the 
pure gauge field A^^^^.^, using operators similar to the operators of specified twist in the opera- 
tor product expansions. These expressions are apparently frame-independent. Based on these 
expressions, we suggest gauge invariant decompositions of the nucleon spin in section 4.1. We 
then pay special attention on the photon spin operator in section 4.2. When restricted on the 
Abelian case, we show that our proposed photon spin operator can describe the spin of light 
modes with a fixed frequency, such as the Laguerre- Gaussian laser modes. We further discuss 
the parton distribution functions constructed from Dirac variables [18-21] and their relations to 
the gauge invariant quark momentums in section 4.3. Finally, we give conclusions in section 5. 
We also have three appendices to give more details of the body of our paper. 



2 Gauge Invariant Angular Momentum Currents As Summations 

For gauge theories, the conserved currents constructed from the Noether theorem are generally 
not gauge invariant. However, we show in this section that gauge invariant currents can be 
constructed from the summations of gauge variant currents. We consider the Lagrangian 

C = -\f;^F^^'' + '-[^l^i^id, - zgA,)^ - (9^^ + tgi^A,)^^^P]. (2.1) 

For the Lorentz invariance of this Lagrangian, the corresponding Noether current can be written 

as 

dC dC - dC 

d{df,A^) d{d^ilj) d{df,^) 

with field variations up to first order of the parameters u^i, as 

SLTAf,{x) = -y'^dfiAf,{x) + u^pA'^{x), = 6x^ = u^^x", (2.3) 
Siri^ix) = -'-ooa^S'^^ijix) - y^d^ij{x), (2.4) 

W(a:) = '-^Pix)u^^S''^ - y^dp^{x). (2.5) 

When Eq. (2.2) is written in a manifest way, it gives the canonical angular momentum current 
as proposed by Jaffe and Manohar [22]. The Lagrangian (2.1) is also invariant under gauge 
transformations. For gauge transformations, the corresponding Noether current is given by 

dC ^ . dC ^ , ^ - dC ,^ . 

d[d^Ay) d[d^i)) d{d^tp) 
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with field variations up to first order of tlie gauge parameter 9 as 



5gtA^{x) = D,e{x) = d,9{x) - ig[A„e{x)], (2.7) 

ScTi^ix) = ig9{x)tp{x), ScTi^ix) = -ig'i/j{x)6{x). (2.8) 

Written manifestly, the current in Eq. (2.6) is 

J^^ = -2TT{Ff"'DJ) - g^lj-f'^eij, (2.9) 

= 2Ti{D^F^"'e) - gip-ff^eilj - 2d^TT{F'"'e). (2.10) 

Here and hereafter we use the normalization Tr(r'^T^) = The first two terms of Eq. (2.10) 

give zero by the equations of motion. The third term satisfies d^dyTi{F^^^9) = 0. So the current 
J'qj, is conserved for any values of the field 0{x). We can suppose a parametrization of the field 

6{x) as 

9{x) = -Uapx'^N^ix). (2.11) 
Then the current J'^j, can be summed with J'j^j, to formulate a new current 

Jf^ = + J^^ = Iw.^M'^"^ (2.12) 

^M«/3 ^ M^aP ^ ^^a^ ^ ^M"/3 + M^^fi ^ (^2.13) 

with identifications of quark parts as 



= le'-^^'i^^pl^i^, (2.14) 



M^7^ = ^[^7'^x"(a^ - igN^)ij -{ai — > /3)] + H.C. (2.15) 

and identifications of gluon parts as 

j^^a(s ^ _2Tr [F'^°(A^ _ Ar/3) _ ^'^/^(A" - AT")] , (2.16) 
M^f = 2Tr [F''^{F''"x^ - F^^x")] (2.17) 
+ 2Tr [x^Ff'^D'^iA'' - iV") - {a < — > /?)] 

Here £quark and £giuon are respectively the corresponding gauge invariant fermion part and gluon 
part of the Lagrangian (2.1), and = — ig[Afj^,-]. The above results are identical with 
our previous ones obtained by using a generalized Lorentz transformation [23]. The cause of 
their consistency is simply that the Lie algebra of the generalized Lorentz transformation used 
in [23] is the direct sum of the algebra of the Lorentz transformation and that of the gauge 
transformation. For A^^ = A^, we obtain the gauge invariant decompositions of Ji [7]; For 
we obtain the gauge invariant decompositions of Chen et al. [8]. The feasibility 
of derivation of gauge invariant currents by the Noether procedure is also discussed in [21, 24] 
through different methods. Their results correspond to special cases of our above discussions 
when A^'^ = v4^ure- Harindranath and Kundu [25], Shore and White [26] ever showed that 
the difference between the decomposition of Ji and that of Jaffe and Manohar is just a total 
divergence by imposing equations of motion. Our above results are also consistent with their 
analysis in the case = AP. 
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3 Novel Expressions for the Pure Gauge Field 
3.1 Abelian Case 

A key step of the proposal of Chen et ah is to express the pure gauge field ^pure terms of 
A'^ in order that no extra freedoms are introduced. For several expressions of the pure gauge 
field, see [10-13]. In this section, we attempt to derive expressions for A^urc through a method 
introduced by Delbourgo and Thompson [27]. For an Abelian gauge field, we consider the pure 
gauge field determined by 

d^[A^{x)-dM^)] = 0^ (3.1) 

which can be rewritten as 

d''dM^)=Tix), ^(x) = 9^A^(x). (3.2) 

A series solution based on the inverse derivative operator has been given in [23, 27, 28]. We try 
to find another kind of solution of Eq. (3.2), we make the ansatz 

oo [f ] 

^(^) = 5Z X]a.,n(x%)"+^x"i ■ ■ -x^^-a^, ■ ■ ■a„,_,ja^a^)"^(x), (3.3) 

s=0 n=0 

where the notation [|] stands for the largest integer which is smaller than |, and Us^n are some 
unknown constants to be determined later on. In the ansatz (3.3), s + 1 counts the number of 
the derivatives appearing in the expansions. The terms in the ansatz (3.3) look very similar to 
the operators of specified twist appearing in the conventional operator product expansions. The 
ansatz (3.3) has the simple property, that is, after being acted by the derivative operator twice, 
it is transformed to be 

oo [f ] 

d'-dM^) = J]^c,,„(x''x,)"x"^ ■ • ■x"--a,, ■ ■ ■a„,_,„(a^a^)"j-(x). (3.4) 

s=0 n=0 

It remains a similar structure to Eq. (3.3), but its coefficients are modified to be 

c^ „^ = {n + l){2d + 4s - 4n)as^n + {s - 2n + 2){s - 2n + l)a^,n_i (3.5) 
+ 4(n + l)as-i,n + 2(s - 2n + l)a^_i,„_i + as-2,n-i, 

with the constraints that ag^n = for n > [|] + 1 ors< —1 orn< —1, which are apparently 
absent from Eq. (3.3). In the above, d is the dimension of the space-time, and hereafter we shall 
work with d = 4. Therefore, a consistent solution of Eq. (3.2) exists if 

co,o = 1, is = 0) (3.6) 
Cs,n = 0. (s>l) (3.7) 

Eq. (3.7) should be satisfied for s > 1. It is a recursive identity, by which ag^n with larger s 
and n can be determined by as^n with smaller s and n. By this recursive identity, the first four 
equations for Qg^n are 



8ao,o 


= 1, 


{s 


= 0,n 


= 0) 


(3.8) 


12ai,o + 4ao,o 


= 0, 


{s 


= l,n 


= 0) 


(3.9) 


16a2,o + 4ai,o 


= 0, 


is 


= 2,n 


= 0) 


(3.10) 


24a2,i + 2a2,o + 2ai,o + ao,o 


= 0. 


is 


= 2,n 


= 1) 


(3.11) 
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Here we have used the relations that ag^n = for n > [|] + 1 or s < — 1 or n < —1. For the 
equations derived from the recursive identity (3.7), the unknown variables match the number of 
equations. So the solution of Eq. (3.1) can be determined by the recursive equations such as 
Eqs. (3.8)-(3.11). The pure gauge field are hence obtained as 

A^^^,,{x) = d^if{x). (3.12) 

From Eq. (3.3), we see that A'^^^^{x) = for the Lorentz gauge d^Afj_{x) = 0. However, the pure 
gauge field determined by Eqs. (3.3) and (3.12) has a complicated formulation, and its transfor- 
mation under gauge transformations is not transparent. In appendix A, we give arguments that 
it has the transformations as as we required. In appendix B.2, we present another solution 
of Eq. (3.1) through a slightly different method, that is, 

i^,,,(x) = A^{x) - ^x«F"^(x) + ^x„,x„,a"iF"^'^(x) (3.13) 

_|_ Ir/^r f) F°'^°'^(t] -\- —tPt f) F^"' ( r] -i- ■ ■ ■ 
8 lb 

By this formulation, At^^^^{x) apparently transforms as A'^{x) under gauge transformations. 
Eq. (3.13) is expected to be equivalent with Eq. (3.12) by similar arguments in appendix A, 
albeit a straightforward proof of their equivalence, which means that we attempt to derive 
Eq. (3.13) from Eq. (3.12) by the method of simple decompositions and combinations, is not 
easy to implement. 

3.2 Non-Abelian Case 

For non-Abelian cases, we need to solve the equation 

Df^iA^ - -U-^d^U) = 0, U{x) = exp-*^<^("), (3.14) 

where (f){x) = (p°-{x)T°- taking values in the Lie algebra of SU{N) and = — ig[A^, ■] stands 
for a covariant derivative. Expanded as a series, the pure gauge field is 

^Cure = -U-'d^'U = d^c^ - '-g[dy, 0] - l9'[[d^<P, + (3.15) 

Because of its non-linearity, Eq. (3.14) is not easy to solve. We attempt to solve it by perturbative 
methods. We suppose Eq. (3.14) has a series expansion as 

= 0(0) + c^0(i) + /0(2) ^ . . . _ ^3_^g^ 

Then Eq. (3.14) can be solved by setting terms of the same power of the coupling g to be zero. 
For the term independent of g, we obtain 

a^a^0(°)(x) = = aM^(x). (3.17) 

This equation has the same structure with Eq. (3.2). So we can solve it by the ansatz (3.3). The 
solution is 

OO [5] 

0(o)(x) = 5^^a,,„(x%)"+V°i ■ ■ ■ ■ ■a„,_,„(a^a^)"^(°), (3.18) 

s=0 n=0 
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with 0'^°) and Aa{x) taking values in the Lie algebra of SU{N), and the coefficients as^n solved 
by Eqs. (3.6) and (3.7). For the term of g of power one in Eq. (3.14), we obtain 

a/^a^<^a) = ^w, (3.19) 

= (3.20) 

In Eq. (3.20), t^'-^^ is known because (/)^°^ has been given by Eq. (3.18). This equation also has 
the same structure with Eq. (3.2). We can solve it by the ansatz (3.3) with replaced by 
The solution is 

oo [f ] 

'/'^'H^) = E E cisA^v^'x'"' ■ ■ ■ ■ ■ ■ d^^^.^id'^d^r^^'K (3.21) 

s=0 n=0 

For the coefficient of in Eq. (3.14), we obtain 

^(2) = S'^ 9^0(2)^ ^322) 

^(2) = 1 [9^0(0), + + l[[a^0W,0{o)],(/>(o)] (3.23) 

2 2 6 

This equation can be solved as above. The solution is 



s=0 n=0 



By the foregoing recursive procedures, a solution of Eq. (3.14) can be derived for the non-Abelian 
case. The pure gauge field Apu^e hence can be obtained for the non-Abelian case. It is apparent 
that ^pui.c(^) = in the Lorentz gauge d^Afj^{x) = by the foregoing constructions. Another 
solution of Eq. (3.14) can be derived by following the procedure in appendix C, which is expected 
to be a generalization of Eq. (3.13) 

A>;^Jx) = A^{x) - ^x,F"^(x) + ^x,,x«,D"^F°^^(x) (3.25) 

"t~ 7^ ai D F (x) + ——X^XpIDq^FP (x) + ■ ■ ■ , 

8 16 

which is also expected to be equivalent with the solution (3.15) obtained above, albeit a formal 
proof is not easy to implement. 

So far, we have derived expressions for the pure gauge field in the Abelian case and the non- 
Abelian case. These expressions shall be used in the next section to construct gauge invariant 
operators of the gluon spin and the photon spin. 



4 Gauge Invariant Decompositions of the Nucleon Spin 

4.1 Generation of Gluon Spin Angular Momentum from Surface Currents 

In section 2, we interpreted the gauge invariant current of angular momentum as summations 
of gauge variant ones. We saw that the gauge invariant decomposition of Ji and that of Chen 
et al. can be accommodated into this framework. In the decomposition of Chen et al. [8], there 
is gauge invariant description of gluon spin. However, in the decomposition of Ji [7], there is no 
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gauge invariant description of gluon spin, while only the total gluon angular momentum current 
is gauge invariant. In this section, we attempt to propose another decomposition of the nucleon 
spin. We begin with Ji's decomposition, that is, the quark parts are 

M!;f = le^'^^'H^.i^, (4.1) 

M^f = ^[^p-f''x''{^^ - %gA^)^l) - (a ^ /?)] + H.C. (4.2) 



+ (r/^"x^ - r/^^a;°)/:q,ark, 



and the gluon part is 



= 2Tr [F^'^{F''''x^ - F'^^x'')] (4.3) 
+ {V^'^x^ - r?'^^a;")£gi,on. 

Apparently there is no natural candidate of the gluon spin operator in the above decomposition. 
In order to introduce descriptions of the gluon spin, following the similar proposal of Belinfante 
for the energy-momentum tensor [29], we construct a gauge invariant surface current 

= ]^dp{x'^J^^''^ - x^^7°), (4.4) 
^.a, ^ _2Tr (F^" A^,^, + + F^^ A^^^^) . (4.5) 

Here the physical field ^phys = A^ — AJ^^^.^ has been defined in section 3. Apparently A^^"'^ is 
also a conserved current, that is, 

d^M^g'^" = 0, (4.6) 

because ^^^^^ is antisymmetrical about its indices /i and p. Adding this current to Af^"^, we 
can obtain a new conserved current 

M^'"/^ = M^g""^ + M''"^. (4.7) 

The gluon parts of this new current are 

M^f = -2Tr (F^- A^,^, + F^^AX. + A^^) , (4.8) 
M^f = 2Tr[F'^,(F'^"x^ - F^^x")] + i:i[x^ F^'" D pA^s - x'^ F^^ D ^A^^^^] 

+ Tr[x^D,F^^^^,^, - x-DpF^^A^^J (4.9) 
+ Tr[x°A^,y,DpF'^^ - x^Al^^^DpF^^] + (r/^-x'^ - r;^'^a;")£gi,on. 

Here = — iglA^, ■] and we have used the relation -Dp^phys = as it is the definition of 
^phys section 3.2. The quark parts of this new current remain the same as Eqs. (4.1) and 
(4.2). The above gluon spin operator (4.8) has been discussed in [12] in the light-cone gauge. 
The second line of Eq. (4.9) is reminiscent of the potential angular momentum in [9]. Of course, 
the surface current (4.4) is not the unique current which can be used to construct descriptions 
of the gluon spin. For other choice, we can consider the current 

^ _2Tr(F^"A^,^, + Ff'^AXs + ^"^^^hys) - 2^g^T{A;^^,[AX^, A^.^J). (4.10) 

This choice induces the gauge invariant Chern-Simons current as descriptions of the gluon spin. 
We shall make more discussions of appropriate descriptions of the gluon spin in the subsequent 
two subsections. 
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4.2 Analysis on Laguerre-Gauss Laser Modes 

In this subsection, when restricted on the Abelian case, we discuss how the operator (4.8) can 
be an appropriate description of the photon spin with a fixed frequency, such as the Laguerre- 
Gauss Laser Modes. As shown by Allen et al. [30, 31], the Laguerre- Gaussian laser modes have a 
well-defined angular momentum, and the angular momentum can be well described by classical 
analysis without the need of quantum field theories. The Laguerre-Gaussian laser modes are 
optical modes with a fixed frequency, and they conform to the source-free Maxwell equations 

Vx^.-f.Vx«-.f. (4.1) 

For E and B with the time dependence e"""^*, the above equations can be rewritten as 

B = X E, E=-VxB. (4.12) 

u u 

In the above descriptions, only gauge invariant field strength are needed. Of course, we can also 
describe the Laguerre-Gaussian modes by using the gauge potential = {A^,A). Then the 
field strength E and B can be expressed by 

E = -VA° - —, B = VxA. (4.13) 

To solve the Maxwell equations, we have the freedom to choose a convenient gauge. We consider 
the Lorentz gauge 

'dt 



^^'A^ = ^ + V ■ A = 0. (4.14) 



In the Lorentz gauge, for optical modes with a fixed frequency, and A can both have the time 
dependence e"*"^*, which are consistent with the Lorentz gauge condition (4.14). So Eq. (4.13) 
can be rewritten as 



E = -VA" + tuA, B = W X A. (4.15) 

Now we begin to consider the descriptions of the photon spin. In the Lorentz gauge d^A^ = 0, 
we know that A'^^^^ = from the discussions in section 3.1. So the photon spin operator (4.8) 
in the Lorentz gauge is simplified to be 

M^f = -{F^'^A^ + F^^'A" + F^^A^"). (4.16) 

We define the photon spin vector as 

S' = ]^e''^ j d^xM'^gi^ (4.17) 

when written manifestly, which is, 

S= f d?x{Ex A + BA^). (4.18) 



Because 5, A^ and A all have the time dependence e ^'^^ as known in electrodynamics, the 
average of the photon spin vector over a time period ^ can be given by 

{S) = ^J d^x{E* X A + Ex A* + B*A^ + BA^*). (4.19) 



Here Z* means the complex conjugate of Z. Using the relation of E and B in Eq. (4.12), we 
obtain 

(S^) = 1 I Sx{E* xA + ^xi'* + -Vx E*A^ - - V X E/^*). (4.20) 
2 J u uj 

After integrating by parts and supposing the surface terms to be zero, we obtain 

{S) = - [ d'xlE* X (1+ -VA°) + Ex{A* - -VA°*)]. (4.21) 
2 J OJ OJ 

Then using the definition of E in Eq. (4.15), finally we obtain 



(S) = -- I d^x{E* X E). (4.22) 



This is consistent with the results in [14, 31, 32]. 

We give some comments here. For an Abelian gauge theory, the conventional Chern-Simons 
current (4.16) appears to be the reduced version of the extended Chern-Simons current (4.10) 
after the Lorentz gauge fixing d^A^{x) = 0. We just show that the Chern-Simons current 
works well as descriptions of the spin of optical modes with a fixed frequency. From the above 
discussions, we may get the impression that the conventional Chern-Simons current (4.16) is 
sufficient for descriptions of the photon spin. However, this impression is not right. Although 
the Chern-Simons current (4.16) works well in classical descriptions as above, it does not work 
as descriptions in quantum theories. As shown by Manohar in the two dimensional Schwinger 
model [33], a forward matrix element of the Chern-Simons current is not gauge invariant. How- 
ever, a gauge invariant extension of the Chern-Simons current (4.10) can yield gauge invariant 
results as shown in [23]. Therefore, it is the gauge invariant extension of the Chern-Simons 
current (4.10) works well both in classical and quantum theories. 

Besides, the foregoing analysis apparently only applies to Abelian gauge theories. Because 
of the nonlinearity of non-Abelian theories, we do not have an appropriate definition of gluon 
modes with a fixed frequency. Therefore, for non-Abelian theories, we still can not determine 
whether it is the operator (4.8) or the operator (4.10) that can be considered as appropriate 
descriptions of the gluon spin by the above analysis, because they degenerate into the same 
formulation when applied to an Abelian theory. 

4.3 Relations to Parton Distribution Functions 

As shown in [7], the quark orbital angular momentum (0AM) can be connected to the gener- 
alized parton distributions. In this subsection, we discuss the relations between the operators 
constructed in section 4.1 and the parton distribution functions. We consider the parton distri- 
bution function 

m) = ^J ^e^'HPSimi^m Xn]i;{Xn)\PS). (4.23) 
Here n'^ = -^(1, 0, 0, 1) is a light-like vector, and U[0; Xn] is a path-ordered gauge link 

U[0; Xn] = Vexp {^9 n^N^\un)du^ , (4.24) 

which makes fq{^) gauge invariant. The quark momentum can be connected to the second 
moment of /g(0 

'dCMO = n,,n,,{PS\mj^H{d^' - z^/Ar^^)^(0)|P5). (4.25) 
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For N^^ = A^, we obtain the gauge invariant quark momentum of Ji [7]; For A^'^ = A'^^^^, we 
obtain the quark momentum considered by Chen et al. [34]. The above definition of parton 
distributions includes the gauge hnk U[0; Xn], which depends on the choice of path (For a dis- 
cussion, see [16]). However, we can use the Dirac variables [18, 19, 21] to formulate another 
definition of parton distributions, which are free of the path dependence. In section 3, in order 
to derive expressions of A^^urc have derived manifest expressions of the unitary matrix U{x), 
which can be used to formulate the Dirac variable 

^^(x) = U{x)ip{x). (4.26) 

By the definition of U{x) in section 3, ipoix) is apparently gauge invariant. The parton distri- 
butions can be constructed from Dirac variables as 

/<f (0 = ^/ ^e^^^{PS\mU~\^)l^U{\n)ii\n)\PS), (4.27) 

which has the bi-local formulation '0£)(a:)'0D(O), and each term of this bi-local expression is 
separately gauge invariant. Its second moment gives 

j daf^ii) = n,,n,,{PS\mY%d^' - ig ■ -U-^d^^U)m\PS). (4.28) 

It shows that the pure field Apuj-c — -U~^d'^U appears in the covariant derivative. So the gauge 
invariant extension of the canonical quark momentum based on the pure gauge field is closely 
connected to the parton distribution based on Dirac variables. From the definition of A'^^^^ in 
section 3, we saw that the quark momentum of Chen et al. can be obtained by subtracting a 
series operator contributions from that definition of Ji. 

When the above applied to the gluon distributions, we have the gluon distributions based 
on Dirac variables 

(0 = ^ / ^e^''HPS\Um^-m-\m{>^n)F^^{\n)U-\\n)\PS). (4.29) 

Here the Dirac variable U{x)F^^{x)U~^{x) is gauge invariant by the definition of U{x) in sec- 
tion 3. The first moment of f^{C,) is 

AG = I dUl^iO = ^{PS\K^\PS). (4.30) 

Here is the Chern-Simons current which satisfies d^K^ = ^F^j^F^^'^ [35]. Because f^iO 
is gauge invariant, we can regard K'^ as the gauge invariant extension of the Chern-Simons 
current. As shown in [23], the gauge invariant extension of the Chern-Simons current also 
satisfies = ^F^^F^*". From this regard, it appears that the gauge invariant extension of 

the Chern-Simons current works as an appropriate description of the gluon spin. 

5 Conclusions 

We have discussed the feasibility of gauge invariant decompositions of the angular momentum 
of the gauge field. We examined the structure of the angular momentum current through the 
Noether theorem. The gauge invariant angular momentum currents are shown to be summa- 
tions of gauge variant currents, which are conserved Noether currents induced separately by the 
Lorentz transformation and the gauge transformation. This summation formulation suggests 



-10- 



that a gauge invariant measurement should include not only the conventional canonical current 
but also the gauge current. We construct novel expressions of the pure gauge field by using a 
tower of operators similar to twist operators in the operator product expansions. These expres- 
sions show that the pure gauge field includes contributions from an infinite operator series. It does 
not seem easy to sum this tower of operators when performing phenomenological calculations. 
We also discussed the appropriate operator description of light modes with a fixed frequency. 
The gauge invariant extension of the Chern-Simons current is shown to yield results consistent 
with other classical analysis. Regarding its frame-independence and its gauge-independence, the 
gauge-invariant Chern-Simons current could be a proper candidate for descriptions of the gluon 
spin, albeit its physical meaning is not transparent because of its manifest non-local formulation. 
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A Transformation of the Pure gauge Field 

In this appendix, we give arguments that the pure gauge field ^pure derived in section 3.1 
transforms as A^^ under gauge transformations. For a gauge transformation, A^^ transforms as 

6A^{x) = d^A{x). (A.l) 

Then by Eq. (3.3), tp transforms as 

OO [f ] 

Sifix) = ^^a,,„(x%)'^+V^---x"--9,,---a„,_,„(9''9;3)"'^-F(x), (A.2) 

s=0 n=0 

SJ^ = 0^5 Af" = d^df'Aix). (A.3) 

From the discussions in section 3.1, it is easy to understand that 6(f in Eq. (A.2) satisfies the 
equation 

d^d^'Sifix) = df,d^A{x). (A.4) 

So Sf can be solved as 

6ip{x) = h{x) + A(x), df^d^hix) = 0. (A.5) 

Here h{x) is independent of A{x). So 6(f{x) equals to A{x) up to a harmonic function. Moreover, 
from Eq. (A.2), we know that 

A{x) = ^ 6ip{x) = 0, (A.6) 

which means that h{x) should be zero and actually 6ip{x) = A{x). So the pure gauge field 
^pure = d^(p{x) transforms as A^. 

B Another expression of the Pure Gauge Field: Abelian Case 

In this appendix, we are going to derive another expression for the pure gauge field, which satisfies 
Eq. (3.1). This expression looks different from Eq. (3.3), but is expected to be equivalent to 
Eq. (3.3). In order to implement this goal, we first solve the pure gauge field through the 
Fock-Schwinger gauge condition, which then can be used to find another solution of Eq. (3.1). 
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B.l Fock-Schwinger Gauge in Abelian Case 

For an Abelian gauge field, we consider a pure gauge field determined by the Fock-Schwinger 
gauge condition 

x^'[A^{x)-d^^Fs{x)] = Q. (B.l) 
To solve Eq. (B.l), we make the ansatz 

oo 

ipFs{x) = 6ix"A,(x) + ^fenx"^ ■ ■ -x^-a,, ■ ■■d^^_,A^Sx)- (B.2) 

n=2 

This expansion is similar to but different from that in section 3.1. A consistent solution of 
Eq. (B.l) exists if 

61 = 1, 62 = -^, 6„ = ?2 = 3,4,5---. (B.3) 
z n 

The solution Eq. (B.2) is not very concise at first sight. However, compared to those expressions 
in section 3.1, Eq. (B.2) has the nice property that, through some simple combinations, the pure 
gauge field that determined by it can be concisely expressed as 

^ure = d'^Fsix) (B.4) 

' ^-^ In + 1)! 



n=2 



where F°'^[x) = d°'A'^{x) — d^A'^lx) is the Abelian field strength, and then the corresponding 
physical field can be given as 

•^phys = A^ — ^pm-j, (B-5) 

2 + 
which satisfies the Fock-Schwinger gauge condition 

^.-K^ys = 0- (B.6) 

Similar but slightly different expressions of the physical field have been given in [36, 37] through 
different methods. From the above, we saw that the pure gauge field (B.4) that determined 
by the Fock-Schwinger gauge condition has a very concise expression and its transformation 
under gauge transformations is also transparent. These concise expressions can be used in next 
subsection to obtain concise expressions of the pure gauge field determined by the Lorentz gauge 
condition. 

B.2 Expression of the Pure Gauge Field in Abelian Case 

Now we try to solve Eq. (3.1) in another way. We decompose <^{x) in Eq. (3.1) as 

ip{x) = (pFs{x) + <p{x). (B.7) 
Then by Eqs. (B.4) and (B.5), Eq. (3.1) can be rewritten as 

d>^dM^) J= = d,M;^^^. (B.8) 



Here ^phys given by Eq. (B.5). This equation has the same structure with Eq. (3.2). It can 
be solved by the ansatz Eq. (3.3), that is, 



^{^) = z2z2 asA^'^pT'''^''' ■ ■ ■ x'-^-'-d^, ■ ■ ■ a„,_,„ {d^d^r^. (B.9) 

s=0 n=0 

In the Abelian case, J-" is apparently gauge invariant by its definition in Eq. (B.8) and by 
Eq. (B.5). So (^(x) is also gauge invariant. By the above solutions, Apure can be written as 

= A^{x) - Ix^F^'^ix) + lx^,x^,d'''F''^''{x) (B.IO) 
2 6 

~l~ T^x^ x^^.Oqi^F (x) -|- x^ X pdctF^ (x) -|- ■ ■ • . 
8 lo 

From the above, we saw that the pure gauge field is decomposed into two parts: w4.[,'^J,p and 
yi{^^j,g. The part ^pu^e determined by the Fock-Schwinger gauge condition, and it has a very 
concise expression (B.4). The part Aj^uj-e is determined by Eq. (B.9), which depends on ^p^ys 
in Eq. (B.5). That is the above decomposition makes the solution (B.IO) to be possible. Here 
we only write out the terms of two derivatives, albeit a expression of all orders is also manifest 
by the above constructions. This new formulation makes the gauge transformation of the pure 
gauge field transparent. 

C Another expression of the Pure Gauge Field: Non- Abelian Case 
C.l Fock-Schwinger Gauge in Non- Abelian Case 

Similar to the discussion of Fock-Schwinger gauge condition in the Abelian case, for non- Abelian 
cases, we need to solve the equation 

x^'iA^ - -Upld^Ups] = 0, Ufs{x) = exp-^^'^^«(^), (C.l) 

where (j){x) = (f)"'{x)T"- taking values in the Lie algebra of SU{N). Expanded as a series, the 
pure gauge field is 

Because of its non-linearity, we attempt to solve Eq. (C.l) perturbatively. We suppose Eq. (C.l) 
has a series solution 

</'F5 = + + 9'<P% + ■■■■ (C.3) 

Then Eq. (C.l) can be solved by setting terms of the same power of the coupling g to be zero. 
For the term independent of g, we obtain 

x^[A,-d,cpP,]=0. (C.4) 

This equation has the same structure with Eq. (B.l). So we can solve it by the ansatz (B.2). 
The solution is 

oo 

(x) = hx^A^ix) + ^ 6„x"^ ■ ■ ■ x'^-a,, ■ ■ ■ 9,„_,A„„(x), (C.5) 



with (j)pl, and Aa{x) taking values in the Lie algebra of SU{N), and the coefficients bn given by 
Eq. (B.3). For the term of g of power one in Eq. (C.l), we obtain 

x'^l^W - d,4l] = 0, (C.6) 

'^[d'^<l>Ps,<l>Ps] = <'- (C.7) 

is known because is given by Eq. (C.5). This equation also has the same structure 
with Eq. (B.l). We can solve it by the ansatz (B.2) with replaced by The solution is 

oo 

(t>%{x) = + ■ ■ -x^-a,, ■ ■ -a^^.^^Hx). (C.8) 

n=2 

For the coefficient of g"^ in Eq. (C.l), we obtain 

x'^K;2)-a,0g)] = o, (C.9) 

'^[d'^Ps^ 41] + lld^Ps^ <l>Ps] + lw<l>Ps. 4l].4l] = (C.IO) 
This equation can be solved as above. The solution is 

oo 

(P%ix) = + ^ . . . ^a„g^^ . . . d^„_^S^l\x). (C.ll) 

n=2 

By the foregoing recursive procedure, a solution of Eq. (C.l) can be derived for non-Abelian 
cases. We saw that the solution of (f)Fs has a very complicated formulation. However, the 

expression for A'^^^^ and ^p^ys ^^^y concise. They are straightforward generalizations of 
Eqs. (B.4) and (B.5), and we have 

^l^ure = '-Up'sd^Ups = A^^ix) - ^x^F^>'{x) + ^x^,x^,D^^F^'>^{x) + ■ ■ ■ , (C.12) 
^^hys = A'- ^Uro = l^^F'^^ix) - lx,,x,,D°^F°^^(x) + ■ • • , (C.13) 

where F°'^^{x) = d°'A^^{x) — d'^A°'{x) — ig\A'^{x),A^[x)\ is the non-Abelian field strength and 
j~)OL _ Q'^ — ig\ji^°i^{^x\ ■]. The physical field here also satisfies the Fock-Schwinger gauge condition. 

C.2 Expression of the Pure Gauge Field in Non-Abelian Case 

In this appendix, we try to solve Eq. (3.14) in another way, similar to the Abelian case, we 
decompose 0(x) as 

0(x) = 0F5(x) + 0(x). (C.14) 

In the series formulation, it is 

0(0) (:c) = cl>fs{x) + 0(°)(x),</>«(x) = + ~4\^). (C.15) 

= <^g)(:,) + 0(2)(^),... . 

Then another solution of Eq. (3.14) can be derived through the similar procedures in ap- 
pendix B.2 and section 3.2. 
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